Bose-Einstein Condensates in Optical Quasicrystal Lattices 
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We analyze the physics of Bose-Einstein condensates confined in 2D quasi-periodic optical lattices, 
which offer an intermediate situation between ordered and disordered systems. First, we analyze the 
time-of-flight interference pattern that reveals quasi-periodic long-range order. Second, we demon- 
strate localization effects associated with quasi-disorder as well as quasiperiodic Bloch oscillations 
associated with the extended nature of the wavefunction of a Bose-Einstein condensate in an optical 
quasicrystal. In addition, we discuss in detail the crossover between diffusive and localized regimes 
when the quasi-periodic potential is switched on, as well as the effects of interactions. 
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The last few years have witnessed a fastly growing in- 
terest on ultracold atomic gases in laser-generated pe- 
riodic potentials (optical lattices, OLs). These present 
neither defects nor phonons, offering a powerful tool for 
investigating the quantum behavior of periodic systems 
under unique control possibilities. Thus, ultracold atoms 
trapped in OLs show fascinating resemblances with solid- 
state physics, which range from Bloch oscillations P, Q 
and Wannier-Stark ladders [3|, to Josephson arrays of 
Bose-Einstein condensates (BECslQ, or to the super- 
fluid to Mott- insulator transition [§J ■ 

These remarkable experiments have been performed in 
regular cubic OLs. However, these lattices do not exhaust 
the rich possibilities offered by optical potentials. More 
sophisticated lattice geometries have been proposed, as 
honey-comb 6] or Kagome and triangular 7] lattices. Be- 
yond, controlled defects may be introduced to generate 
random or pseudorandom potentials [j| , allowing for the 
realization of Kondo-like physics |9j, Anderson localiza- 
tion and Bose-Glass phases |T(| . Exploiting this possibil- 
ity, laser speckle fields have been emplo yed ver y re cently 
to produce BECs in random potentials [III E2> 03 open- 
ing very exciting new experimental possibilities. 

Bridging between ordered and disordered structures, 
quasicrystals (QC) have attracted a wide interest since 
their discovery in 1984 01 . QCs are a long-range ordered 
materials but without translational invariance, and con- 
sequently they share properties with both ordered crys- 
tals and amorphous solids |15| . In particular, QCs show 
intriguing structure |l(| as well as electronic conduction 
properties 0] at the border between ordered and disor- 
dered systems. 

Surprisingly, up to now, only few works have been de- 
voted to optical analogues of QCs, despite of the fact 
that OLs offer dramatic possibilities for designing a wide 
range of geometries [l8|. Optical QCs have been first 
studied in laser cooling experiments ^||, in which the 
atomic gas, far from quantum degeneracy, was confined 
in a dissipative OL where quantum coherence was lost 
due to spontaneous emission. In these systems, the tem- 
perature and spatial diffusion were found to behave simi- 



larly as in periodic OLs. The physics of ID quasiperiodic 
OLs has also been subject of recent research in the con- 
text of cold atomic gases, including a pro pos al for the 
atom-optical realization of Harper's model [2(j, and the 
analysis of Fibonacci potentials |2lJ . 

In this Letter, we study the dynamics of a BEC in a 
2D optical QC. First, we show that the BEC wavefunc- 
tion displays quasi-periodic long-range order, a property 
that may be easily probed via matter-wave interferome- 
try. Second, we show that macroscopic quantum coher- 
ence dramatically modifies the transport on the lattice 
compared with the dissipative case. On the one hand, 
due to quasi-disorder in optical QC, spatial localization 
occurs, in contrast to ballistic expansion in periodic lat- 
tices. The crossover between ballistic expansion and lo- 
calization is analyzed when the quasi-periodicity of the 
lattices is continuously increased. On the other hand, we 
show that due to extended character of the BEC wave- 
function, Bloch oscillations take place. These oscillations 
are however quasi-periodic rather than periodic. Addi- 
tionally, we briefly discuss the effects of the interatomic 
interactions in the BEC diffusion. 
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FIG. 1: Left: laser arrangement (see text for details). Right: 
Quasi-periodic lattice potential for Nb — 5, Vo < and ay = 
0. White points correspond to potential minima. 



In the following we consider a dilute Bose gas trapped 
in the combination of a smooth harmonic potential 
VhoC?) = f (wi^i+wf* 2 ) plus an OL Matt(^±). 
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In the previous expression, M is the atomic mass, u>j are 
the harmonic trap frequencies, and ~r* ± = (x, y) is the 
position vector on the lattice plane. We assume lo z to be 
large enough to keep a 2D ph ysic s on the xy plane. We 
consider a laser configuration |l9j consisting on N\, laser 
beams arranged on the xy plane with TVb-fold symmetry 
rotation (Fig. The polarization ~e j of laser j with 

wavevector k j is linear and makes an angle ctj with the 
xy plane. The optical potential is thus |22| : 



Matt(^±) = 



V 



e j e 
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where < £j < 1 stand for eventually different laser 
intensities and ipj are the corresponding phases. In the 
following, we are mostly interested in the fivefold sym- 
metric configuration (iVb = 5, £j = 1), similar to the 
Penrose tiling 23], which supports no translational in- 
variance (see Fig.^l. The lattice displays potential wells 
which are clearly not periodically arranged. We also con- 
sider the configuration obtained by switching off lasers 1 
and 4, which results in an anisotropic periodic lattice. 

1 ) Equilibrium properties - The stationary BEC 
wavefunction ipo is obtained from evolution in imaginary 
time of the 2D Gross-Pitaevskii equation (GPE): 
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where <?2D = \ZSTrh 3 uj z /Ma sc , with a sc the s-wave scat- 
tering length. In order to elucidate the long-range order 
properties of the BEC, we compute the momentum distri- 
bution of ip . In the periodic case (Fig.|!Ji), as expected, 
the momentum distribution displays discrete peaks cor- 
responding to combinations of elementary basis vectors 
of the reciprocal lattice, nilti + ri2~K2 with integer co- 
efficients n i and ri2- As obtained in previous experimen- 
tal works |24l |25| , this confirms the periodic long-range 
order of ipo. The quasi-periodic case (Fig. is more 
intriguing, resulting in a more complex structure. The 
momentum distribution also displays sharp peaks, being 
the signature of a long-range order which is quasiperi- 
odic rather than periodic |26| . As in the periodic case, 
the positions of the peaks are linear combinations of inte- 
ger numbers of Nb = 5 wavevectors: Ylj=o 1 n i~^i wnere 

~kq = k\— ko and ~~k j — lZ[(f>j]~Ko is the wavevector 
obtained by a rotation of angle <pj — 2nj/Nb of ~kV The 
reciprocal lattice thus clearly shows a five-fold rotation 
symmetry incompatible with any translation invariancc 
[2?| . This resembles the Penrose tiling |2^| and the first 
solide state QCs observed via Bragg diffraction [Tij . 

The discussed momentum distribution can be directly 
imaged via matter-wave interferometry after a time-of- 
flight expansion [28j. Indeed, although the interactions 
are crucial for determining local populations of each po- 
tential well, they do not contribute significantly to the 
free BEC expansion after release from the trap [24( . Such 



measurements, standard in periodic OLs [24j,|25(, can be 
easily extended to quasi-periodic ones. 
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FIG. 2: Matterwave interference pattern of a BEC released 
from a combined OL and harmonic trap: a) periodic case; 
b) quasi-periodic case. Both correspond to 87 Rb and Vo = 
— 10-Br, where En = h 2 k 2 /2M is the recoil energy. 

2) Quantum transport - Certainly, not all physical 
properties of optical QCs can be directly interpolated 
from the behavior of periodic lattices. Indeed, solid QCs 
show intriguing dynamical properties that are not yet 
completely understood 15]. In the following, we investi- 
gate dynamical properties of quasi-periodic lattices. 

Coherent diffusion - Starting from the equilibrium 
wavefunction ipQ: we consider the situation in which the 
harmonic trap is switched-off at t = 0, letting the BEC 
evolve in the OL. The BEC expansion is then computed 
using a Crank-Nicholson algorithm for the real time- 
dependent GPE J5J. Figure |3 (inset) shows the time 
evolution of (x 2 ) and (y 2 ) of the interacting BEC along 
x and y respectively. In the periodic case, the conden- 
sate expands coherently as one expects from tunnel cou- 
plings between adjacent lattice sites. In order to infer 
a convenient fitting functional for the expansion of the 
interacting BEC, let us recall that in free space for large 
times {r 2 (t)) oc v 2 t 2 with Vj oc l/M |2^. In periodic 
lattices, the inertia is enhanced, and the expansion is ex- 
pected to be as in free space but substituting the atomic 
mass M by an effective mass M* > M . We thus expect 
Vj oc l/M* and M/M* oc J where J is the site-to-site 
tunneling rate pflj . Numerical computations for various 
depths of the lattice potential Vq show that v 2 decreases 
exponentially with Vq as expected from the well-known 
exponential decay of J. Anisotropic ballistic expansion 
of the BEC reflects the anisotropy of our periodic lattice. 

The behavior of the BEC in the quasiperiodic lattice 
is dramatically different, since after a short transient 
the BEC localizes ^ (inset of Fig. 0). This behavior 
strongly contrasts with results obtained in the context 
of laser cooling where similar classical normal expansion 
((r 2 (t)} ~ 2D 7 t) was found for both periodic and quasi- 
periodic OLs [13 • Here, spatial localization is a coherent 
effect induced by quasi-disorder due to the lack of period- 
icity. Indeed, the BEC populates localized (Wannicr-like) 
states centered on each lattice site. In the periodic case 
these states have all the same energy and are strongly 
coupled through quantum tunneling. On the contrary, 
in the quasi-periodic lattice, the sites have different en- 
ergies. In particular, the typical difference of depths of 
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adjacent sites (denoted A below) can be of the order of 
magnitude of (but smaller than) the potential depth. The 
tunneling is not resonant and the BEC localizes. 

The remarkable flexibility of OLs 01 allows for the ac- 
curate study of the competition between tunneling and 
quasi-disorder. By ramping up gradually the intensity 
of lasers 1 and 4 while keeping constant 0, 2 and 3, one 
turns continuously from an anisotropic periodic lattice to 
a five-fold symmetric quasiperiodic one, and hence from 
ballistic expansion to spatial localization. For small in- 
tensity of the control lasers 1 and 4, the quasiperiodicity 
is mainly compositional (the sites are still periodically 
displaced but the on-site energies are different from site 
to site) . We define the quasi-disorder A parameter as the 
variance of the differences of on-site energies in adjacent 
sites. From the previous discussion and to compare to 
the results of the non-degenerate case [l^ we fit 



(r 2 J (t)) = (r%) + 2D ] t + vh 2 



(3) 



In the considered range of parameters, all calculations 
fit well with Eq. @ with a negligible diffusive term 
2Djt. We characterize the expansion along the in- 
direction through the ballistic velocity Vj. The behav- 
ior of v x versus the quasi-disorder parameter A is shown 
in Fig. for the interacting BEC and it is compared 
to the non-interacting case. For the latter, we simulta- 
neously switched off the interactions at t = In 
both cases, as expected, coherent diffusion dramatically 
decreases when quasi-disorder increases. Spatial local- 
ization occurs for A > J Xl with J x the tunneling rate 
between adjacent sites along the x-direction 33] . This 
supports the interpretation that competition between co- 
herent tunneling and inhomogeneities turns into localiza- 
tion as soon as tunneling becomes non-resonant. 
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FIG. 3: Crossover from ballistic to localization regimes with 
(squares) and without (circles) interactions, for Vb = — 7.5i?R. 
Inset: Coherent diffusion in periodic and quasi-periodic lat 
tices for Vo = -5£r. Fits to (r 2 (t)) 
the periodic case are also shown. 
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To understand the effect of interactions that can help 
(A < 0.2J X in FigGU or hinder (A > 0.2J X in FigG§ 
diffusion, note first that two phenomena contribute to 
localization: (i) initial inhomogeneities (due to disorder 



and harmonic confinement) that appear in the dynamics 
through the interaction term g2D\4>o\ 2 [see Eq.@] and 
(ii) inhomogeneities associated to quasi-disorder. Be- 
cause of these inhomogeneities, quantum tunneling is 
not resonant and thus less efficient. However, during 
diffusion, the interaction energy is converted into ki- 
netic energy and this tends to fasten the expansion. For 
small quasi-disorder, the second phenomenon dominates 
so that interactions contribute to expansion whereas 
for larger quasi-disorder, the inhomogenities significantly 
hinder tunneling so that interactions contribute to local- 
ization. The non-trivial interplay between disorder and 
interactions will be the subject of further research. 

Quasiperiodic Block oscillations - One of the most 
appealing predictions of the quantum theory of solids 
|27| is that homogeneous static forces induce oscilla tory 
rather than constant motion in periodic structures |34| . 
The corresponding Bloch oscillations have already been 
observed in superlattice superconductors [35j and on cold 
atoms in OLs 0,0- It is a fundamental question whether 
such a phenomenon also exists in less ordered systems like 
QCs. Arguments based on general spectral properties of 
QCs |3a and numerical simulations of ID Fibonacci lat- 
tices [37J support the existence of Bloch oscillations in 
quasiperiodic lattices. However, to the best of our knowl- 
edge, this effect has never been observed experimentally. 
Using accelerated lattices 0, or gravity [38| (we con- 
sider the latter), this question can be addressed experi- 
mentally in the discussed arrangement (Fig. Starting 
from ipo we switched-off the harmonic trap and the inter- 
actions at time t = and tilt the quasiperiodic lattice in 
the x direction 39] . The latter evolution of the quantum 
gas is shown in the inset of Fig. 0] We find noisy-like 
oscillatory motion in the (tilted) x direction and no mo- 
tion in the (non-tilted) y direction. The oscillations in 
the x direction are clearly not periodic 40]. However, 
they definitely have an ordered structure, which is evi- 
denced by the appearance of discrete sharp peaks in the 
time Fourier transform of the BEC mean position (x( t)) 
(Fig. corresponding to a quasiperiodic motion [26j. 

The Bloch-like quasiperiodic oscillations can be inter- 
preted as follows. Both in periodic and quasiperiodic 
lattices, the BEC wavefunction extends over many lat- 
tice wells, and can be decomposed into a sum of localized 
(Wannier-like) states. Due to the applied external force, 
these energy states are arranged in a Wannier-Stark lad- 
der. In periodic lattices, the energy separation AE be- 
tween the ladder states is fixed, leading to periodic Bloch 
oscillations of period oc h/ AE. However, for quasiperi- 
odic lattices, a discrete set of different (non commen- 
surate) differences of on-site energies in adjacent wells 
occurs (i.e. a non-equally spaced Wannier-Stark ladder) 
leading to quasiperiodic (instead of periodic) oscillations. 
Purely random potentials would result in a continuous set 
of differences of on-site energies leading, as expected, to 
the disapearance of any sort of Bloch-oscillations. 

Summarizing, we have investigated the physics of 
BECs trapped in optical QCs. We have shown that 
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FIG. 4: Fourier transform of the mean position of a BEC in 
a periodic or quasiperiodic lattice. Inset: Time evolution of 
a BEC in a tilted quasiperiodic lattice. All beams have the 
same intensity, Vb = — 2_Er and Vtiit = 0.002_Br x kx. 

(i) the equilibrium BEC wavefunction displays long- 
range quasiperiodic order and that (ii) quantum trans- 



port shares properties with both ordered and disordered 
systems. On the one hand, because of quasi-disordered 
inhomogeneities, diffusion turns from ballistic to localiza- 
tion when quasi-periodicity is switched on. On the other 
hand, because of coherence extending over several lattice 
sites, quasiperiodic Bloch oscillations occur in quasiperi- 
odic BECs. 

The discussed arrangement can be easily generated 
using standard techniques, offering an exciting tool 
for controlled studies of the transition between peri- 
odic, quasiperiodic and fully disordered systems in cold 
gas es, a major topic of current experimental research 
{Til 111. IT^j. In addition, the system can be used to ad- 
dress experimentally some unsolved issues on QCs, as e.g. 
the application of the renormalization theory to the 2D 
case |l5j . 
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